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INTRODUCTION 


This  project  is  a  continuation  of  the  earlier  ARO  Grant  DAAG-29-79-G-0017. 
The  overall  objective  of  this  project  was  to  obtain  additional  experimental  data  on  peri¬ 
odic  turbulent  boundary  layers  subjected  to  zero,  as  well  as,  adverse  pressure  gradients, 
from  the  unsteady-flow  water  tunnel  constructed  under  the  earlier  ARO  support.  These 
experimental  data  would  not  only  provide  a  comprehensive  data  base  for  use  by  the 
research  community  but  also  would  directly  be  useful  in  the  development  of  a  theoret¬ 
ical  framework  for  the  description  of  unsteady  turbulent  boundary  layers.  The  present 
contract  work  was  started  at  the  University  of  Iowa  (UI)  on  October  1,  1982.  While 
the  work  was  in  progress,  the  PI,  Dr.  B.R.  Ramaprian  resigned  his  position  at  the  UI 
to  take  up  the  position  of  Professor  of  Mechanical  Engineering  at  Washington  State 
University  (WSU)  with  effect  from  August  16,  1985.  With  the  approval  from  ARO,  the 
research  work  was  continued  by  him  at  WSU,  under  a  subcontract  from  UI  to  WSU, 
with  Professor  F.  Stern  of  the  Department  of  Mechanical  Engineering  at  UI  acting  as 
a  co-principal  Investigator  and  being  responsible  for  reporting  to  ARO.  The  experi¬ 
ments  were  completed  at  Iowa  and  the  unfinished  data  processing  and  analysis  were 
completed  at  WSU. 


Problems  Studied 

The  following  tasks  were  performed  under  this  contract. 

1.  Experiments  on  periodic  turbulent  boundary  layers  in  time-mean  zero  pressure 
gradient. 

2.  Theoretical  Analysis  of  unsteady  turbulent  boundary  layers. 

3.  Computational  studies  on  unsteady  turbulent  boundary  layers. 

4.  Theoretical  analysis  of  the  performance  of  a  surface-mounted  heat-flux  gage  as 
a  skin-friction  gage  in  unsteady  flow. 

5.  Experiments  on  periodic  turbulent  boundary  layers  in  time-mean  adverse  pres¬ 
sure  gradient. 

The  work  performed  and  the  important  results  obtained  are  summarized  in  the 
following  sections.  More  details  can  be  found  from  other  sources  such  as  thesis,  journal 
articles,  interim  reports,  etc.  cited  in  this  report. 


EXPERIMENTS  IN  ZERO  PRESSURE  GRADIENT 


These  experiments  were  performed  on  a  turbulent  boundary  layer  which  was 
subjected  to  a  freestream- velocity  variation  of  the  form 

U,[t)  =  U0{  1  +  c  sinwt)  (1) 

with  the  time  mean  velocity  U0  and  the  relative  amplitude  of  oscillation  e  being  both 
constant. 


Apparatus  and  Instrumentation 


The  unsteady-flow  water  tunnel  built  under  the  earlier  ARO  support  was  used 
for  these  experiments.  Figure  1  shows  the  layout  of  the  tunnel.  It  works  under  a 
constant  head  of  8  m.  The  test  section  is  2.4  m  long,  with  a  rectangular  cross-section 
50  cm  in  width  and  22.5  cm  iu  height.  Its  bottom  wall  used  as  the  test  surface  is  a 
smooth  brass  plate.  Transition  was  promoted  by  means  of  a  strip  of  14-mesh  brass 
screen  (15  cm  wide)  glued  to  the  test  surface  at  a  distance  of  about  22.5  cm  from  the 
end  of  the  contraction.  The  side  walls  of  the  tunnel  are  9.5  cm  thick  steel  plates,  with 
5  windows  whose  centerlines  are  located  at  approximately,  x=48,  69,  90,  142,  and  203 
cm.  These  locations  are  designated  as  stations  1,2, 3, 4  and  5  respectively.  The  top 
wall  of  the  test  section  is  a  Plexiglas  sheet  37.5  mm  thick.  In  addition,  a  6.4  mm 
thick  flexible  Plexiglas  sheet  was  used  as  a  ’false’  wall.  By  adjusting  its  position  to 
correct  for  the  displacement  thickness  of  the  boundary  layer  along  the  walls,  a  nearly 
zero  (time-mean)  longitudinal  pressure  gradient  was  obtained  in  the  test  section.  The 
tunnel  ends  in  a  steel  cylinder  (40  cm  diameter  x  60  cm  long)  whose  downstream  end 
is  closed.  The  cylinder  has  two  longitudinal  rectangular  slots,  30  cm  x  2.54  cm,  located 
180  degrees  apart,  from  which  the  water  exits.  The  area  of  opening  of  the  slots  is 
varied  by  a  rotating  profiled  sleeve  driven  by  a  3  H.P.  geared  D.C.  motor  whose  speed 
can  be  regulated  to  within  1/4%.  Each  complete  rotation  of  the  sleeve  corresponds  to 
two  oscillation  cycles.  The  sleeve  profile  is  contoured  to  produce  a  sinusoidally  varying 
velocity  in  the  test  section  at  the  desired  frequency.  The  details  of  design  of  the  sleeve 
are  given  in  Ramaprian  and  Tu  (1982). 

A  frequency-shifted,  tracker  based,  two-component  Laser  Doppler  Anemometer 
(LDA)  was  used  to  measure  the  instantaneous  longitudinal  and  normal  velocity  com¬ 
ponents  in  the  boundary  layer  of  the  test  surface.  A  surface-mounted  TSI  heat-flux 
gage  (HFG)  operated  by  a  DISA  constant-temperature  anemometer  was  used  for  the 
measurement  of  wall  shear  stress.  A  theory  for  the  use  of  this  probe  in  unsteady  flow 
is  described  in  section  4.  Measurement  of  wall  shear  stress  was  made  at  four  stations: 
x=48  cm,  80  cm,  142  cm  and  203  cm. 

A  HP/1000  minicomputer  was  used  for  data  acquisition  and  processing.  Signals 
from  the  LDA  and  HFG  were  sampled  at  100  fixed  phase  positions  in  the  oscillating 


cycle.  The  instantaneous  data  obtained  were  processed  using  the  well  known  "triple 
decomposition”  principle  (see  for  eg.,  Hussain  and  Reynolds  1970).  This  means  that 
any  instantaneous  property  9  (9  =  U,V  or  r„)  is  expressed  as 

9  =  9{x,y)  +  9p{x,y,6p)  +  <(>(x,y,t)  (2) 

where  9  is  the  time-mean  value,  9P  the  periodic  deterministic  component,  <f>  the  tur¬ 
bulent  fluctuation  and  0P  the  phase  position  within  a  cycle.  The  first  step,  in  practice, 
is  to  obtain  the  ensemble-averaged  value  ( 9)(x,y,0p )  by  averaging  the  instantaneous 
values  at  identical  $p  positions  over  a  large  number  of  oscillation  cycles.  The  time-mean 
value  9  is  then  obtained  by  averaging  (9)  over  the  complete  cycle  (0  <  9P  <  2*).  The 
periodic  component  9P  is  the  difference  (9)  —  9.  Finally,  the  turbulent  fluctuation  <f>  is 
simply  9  —  (9).  This  is  a  statistical  quantity,  for  which  the  phase-averaged  and  time- 
averaged  mean  square  values,  (4>2)  and  <f>*t  respectively,  are  also  calculated  using  the 
above  procedure.  The  instantaneous  velocities  U  and  V  were  ensemble  averaged  in  this 
manner  over  1000  effective  cycles  to  obtain  ( U ),  (V),  (ua),  (v2)  and  (uv)  as  functions  of 
the  phase  position  9P.  For  more  details  of  the  whole  procedure,  see  Ramaprian  and  Tu 
(1982). 


Experimental  Conditions 

The  nominal  mean  velocity  Ua  and  the  amplitude  of  oscillation  t  of  the  free 
stream  were  90  cm/sec  and  0.4  respectively.  Two  different  frequencies,  /  =  u//2x  =  0.5 
Hz  and  2.0  Hz  were  studied.  The  turbulent  bursting  frequency  /fc  in  the  boundary 
layer,  which  can  be  regarded  as  a  characteristic  frequency  of  turbulence  was  estimated 
from  the  criterion  of  Rao,  Narasimha  and  Badri  Narayanan  (1971),  to  vary  from  about 
13  to  4.5  Hz  from  the  first  to  the  last  measuring  station.  It  is  seen  that  the  higher  of 
the  two  oscillation  frequencies  is  comparable  to  the  bursting  frequency,  at  least  for  the 
downstream-most  station. 

Preliminary  surveys  showed  that  U,  varied  by  less  than  5%  along  the  entire  test 
section  length  and  less  than  2%  in  the  range  90  cm  <  x  <  210  cm.  The  amplitude 
e  was  also  found  to  remain  constant  within  the  same  limits.  The  freestream  velocity 
modulation  was  sinusoidal  everywhere  along  the  test  section  with  less  than  1%  total 
harmonic  distortion,  as  confirmed  by  a  Fourier  analysis  of  the  velocity  signal.  It  was 
found  from  several  spanwise  traverses,  that  the  boundary  layer  was  acceptably  two- 
dimensional  over  the  central  30%  of  the  span  (corresponding  to  about  46),  with  respect 
to  both  mean  and  oscillatory  properties.  For  example,  these  properties  were  uniform  to 
within  2%  in  mean  value  and  amplitude,  and  within  0.5  degree  in  phase  in  this  region. 
All  these  initial  experiments  are  documented  in  detail  in  Menendez  and  Ramaprian 
(1983). 


Summary  of  Experimental  Results 


All  the  experimental  data  have  been  stored  on  magnetic  tape  in  the  format 
suggested  by  Dr.  Lawrence  Carr  of  NASA,  Ames  Research  Center,  for  archival  as 
AGARD  data  base.  These  data  can  be  obtained  by  writing  to  either  Dr.  Carr  or  the 
author.  Only  some  typical  data  are  presented  and  discussed  in  this  report. 


Time-Mean  Flow  Properties 

Figure  2  shows  the  variation  of  the  time-mean_yalues  of  Re,  —  Uo0/u,  H  =  6*/ 9 
and  Ct  —  2f m I pUo  with  Re,  =  U0x/u,  where  6*  and  9  are  the  time-mean  displacement 
and  momentum  thickness  respectively.  It  is  seen  that  the  values  of  Re,  and  H  are  the 
same  for  the  two  unsteady  flows  but  are  significantly  different  from  those  for  steady  flow 
at  the  upstream  stations.  The  time-mean  wall  shear  stress,  on  the  other  hand,  is  only 
slightly  affected  at  all  stations  by  unsteadiness.  In  fact,  the  differences  in  Cj  between 
the  three  flows  are  too  small  to  be  clearly  distinguished  from  data  scatter.  The  values  of 
H  and  the  slope  dRe,/dRea  also  indicate  that  the  unsteady  and  steady  flows  approach 
each  other  as  Re,  increases. The  data  trend  suggests  that  the  initial  differences  between 
unsteady  and  steady  flows  are  the  result  of  shifts  in  the  virtual  origin.  In  fact,  a  careful 
study  of  the  data  indicated  (see  Menendez  and  Ramaprian  1083  for  details)  that  the 
steady  flow  is  not  fully  developed  up  to  station  3,  (possibly  due  to  the  ineffectiveness  of 
the  boundary-layer  trip).  It  appears  that  the  imposition  of  oscillation  on  the  flow  has 
the  effect  of  promoting  a  faster  evolution  of  the  flow  (especially  the  outer  layer),  into  a 
a  "standard”  boundary  layer.  The  evolutions  of  the  mean  velocity  distributions  in  the 
usual  inner  and  outer  coordinates  are  shown  in  figures  3(a)-(d).  The  scaling  velocity 
u.  used  in  these  figures  is  obtained  from  the  measured  wall  shear  stress  r„  as  rv/pl/ 3 
(usL’.0  the  fitted  lines  shown  in  figure  2b).  The  scaling  length  6  is  obtained  from  a 
3-parameter  fit  of  the  measured  mean  velocity  distribution  in  the  region  y+  >  40  to 
the  Coles’  "log+wake”  distribution 

V  =  +  5.5  +  _  2,*)]  (3) 

where  k  —  0.41,  rj  —  y/6  and  C/e,ij,6 ,  are  the  three  free  parameters.  It  is  seen  that 
the  velocity  distributions  in  the  inner,  as  well  as  outer  coordinates  approach  their 
respective  steady-flow  distributions  towards  the  downstream  locations.  Once  again,  the 
departure  observed  in  the  upstream  stations  are  the  result  of  the  different  "histories” 
of  the  three  flows  and,  hence,  are  specific  to  the  apparatus.  Some  difference  is  still 
observed  between  the  steady  and  unsteady  flows  at  the  most  downstream  station,  at 
the  oscillation  frequency  of  0.5  Hz  but  not  at  2  Hz.  This  is,  however,  small  and  in  view 
of  the  uncertainties  in  flow  development,  it  is  difficult  to  draw  strong  conclusions  about 
the  effect  (if  any)  of  oscillation  on  the  time-mean  flow. 


The  time-averaged  turbulent  stresses  u3,  v3  and  uv  are  shown  in  figure  4  for  the 
downstream-most  station  (station  5).  The  distributions  at  2  Hz  seem  to  be  very  nearly 
the  same  as  in  steady  flow.  The  distributions  at  0.5  Hz,  however,  show  some  departure 
from  steady  flow.  In  fact,  a  study  of  the  data  from  all  the  stations  indicated  some 
correlation  between  the  effect  on  these  turbulence  properties  and  the  reduced  frequency 
wz/i/,  with  the  effect  being  zero  at  very  small  and  very  luge  values  of  ujx/u  and  being  a 
maximum  around  ux/v  =  5  to  8.  The  authors  have  no  satisfactory  explanation  for  this 
observed  trend.  Other  investigators  have  reported  no  effect  of  imposed  unsteadiness 
on  the  time-averaged  turbulence  properties. 


Oscillatory  Motion 

The  response  of  the  oscillatory  velocity  to  to  the  imposed  oscillation  was  ev¬ 
erywhere  found  to  be  harmonic.  The  total  harmonic  distortion  was  less  than  1  %.  It 
is  therefore  adequate  to  consider  only  the  fundamental  component  of  the  oscillatory 
velocity.  The  response  of  the  other  properties,  however,  may  be  significantly  distorted. 
This  is  seen  typically  from  figure  5  which  shows  the  variation  of  the  ensemble  aver¬ 
aged  wall  shear  stress  (r„)  and  the  Reynolds  shear  stress  (uv)  during  the  oscillation 
cycle  (the  latter  at  three  different  points  in  the  boundary  layer)  at  station  5.  However, 
Fourier  analysis  indicated  that  even  in  this  case,  the  first  harmonic  is  still  the  dominant 
one.  Hence,  only  this  component  will  be  considered  for  discussion  in  this  report.  Com¬ 
plete  information  on  all  properties  (including  all  the  significant  Fourier  components) 
has,  been  stored  on  magnetic  tape.  In  spite  of  the  harmonic  distortion  in  the  (uv) 
variations,  it  is  very  clearly  seen  that  (uv)  is  modulated  very  strongly  near  the  wall 
(and  so  also  is  ( Ct )).  The  modulation  in  (uv)  becomes  weaker  with  increasing  distance 
from  the  wall  until  it  becomes  zero  in  the  outer  part  of  the  boundary  layer  resulting  in 
a  ” frozen”  turbulence  structure  in  that  region.  This  observation  is  in  agreement  with 
other  experiments.  Figure  5  also  shows  that  there  is  a  strong  phase  variation  across 
the  region  where  the  amplitude  modulation  is  significant.  This  region  will  henceforth 
be  referred  to  as  the  "unsteady  vortical  layer”. 

Some  typical  results  for  the  phase  and  amplitude  distributions  of  (17),  (u3)  and 
(uv)  are  shown  in  figures  6  and  7.  These  data  correspond  to  station  5.  The  amplitude 
of  ( U )  exhibits  the  well  known  trend  with  a  slight  overshoot  (7-8  %  of  the  freestream 
amplitude).  The  extent  of  the  region  of  amplitude  variation,  i.e.,  the  thickness  of 
the  unsteady  vortical  layer  decreases  with  frequency.  Beyond  this  layer  is  a  layer  of 
constant  amplitude  or  the  so-called  "slug-flow”  layer.  A  quantitative  measure  of  the 
thickness  of  the  unsteady  vortical  layer  is  the  distance  ym  from  the  wall  to  the  point  of 
maximum  amplitude.  The  quantity  wym3/2i/1/J  which  is  the  ratio  of  ym  to  the  thickness 
of  the  laminar  Stokes  layer  (or  the  unsteady  viscous  layer)  has  a  value  of  about  10,  at 
both  frequencies  studied,  indicating  that  the  unsteady  effects  to  much  larger  distances 
in  turbulent  flow  than  in  laminar  flow.  This  result  is  in  agreement  with  the  pipe- 
flow  data  of  Tu  and  Ramaprian  (1983).  The  phase  lead  of  (C7),  shown  in  figure  6 
(b)  also  exhibits  the  well  known  trend  of  increasing  lead  towards  the  wall.  However, 
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since  data  could  not  be  obtained  closer  than  1  mm  from  the  wall,  the  detailed  manner 
in  which  the  phase  lead  of  the  velocity  approaches  the  phase  lead  of  the  wall  shear 
stress  cannot  be  understood  from  these  measurements.  This  is  not  a  serious  drawback, 
however,  since  extensive  near-wall  data  of  high  quality  are  available  from  the  work  of 
Jayaraman,  Parikh  and  Reynolds  (1982).  The  emphasis  on  the  present  experiments 
was  in  obtaining  the  Reynolds  shear  stress  data. 

The  amplitude  of  oscillation  of  the  turbulence  properties  (u3)  and  (— uv),  nor¬ 
malized  by  2eUj,  is  presented  in  figure  7.  It  can  be  shown,  by  considering  small 
perturbations, that,  with  this  nondimensionalization,  the  resulting  quantities  should  be 
very  nearly  equal,  in  quasisteady  flow  at  small  t,  to  the  steady-flow  values  of  u*/U2 
respectively  (Shemer,  Wygnanski  and  Kit  1985).  Hence,  the  significant  departure  of 
the  distributions  in  figure  7  from  the  steady-flow  distributions  of  figure  4  suggests  that 
the  quantities  (u3)  and  (uv)  do  not  respond  in  a  quasisteady  manner,  at  both  the  os¬ 
cillation  frequencies.  Further,  it  is  seen  that  the  region  of  significant  variation  in  the 
amplitudes  of  these  properties  corresponds  to  the  unsteady  vortical  layer  observed  in 
the  case  of  (17).  Beyond  this  region,  the  amplitude  is  essentially  zero,  indicating  a  state 
of  frozen  turbulence,  in  agreement  with  previous  investigations. 

The  phase  distributions  of  (u3)  and  (uv),  relative  to  the  freestream  velocity,  are 
presented  in  figures  7(c)  and  (d).  These  indicate  phase  lag  in  both  cases,  which  tend 
to  zero  at  or  close  to  the  wall  (where  the  turbulence  is  produced),  and  increase  to  360 
degrees  at  the  outer  edge  of  the  unsteady  vortical  layer  (indicated  by  a  vertical  line, 
for  /  =  2  Hz).  The  data  scatter  observed  beyond  this  layer,  for  /  =  2  Hz,  is  a  result  of 
the  Fourier  analysis  of  the  very  small  (spurious)  modulation  of  velocity  in  that  region 
and  is  not  significant. 

The  amplitude  and  phase  distributions  of  the  turbulence  properties  are  in  gen¬ 
eral  agreement  with  those  measured  by  Ramaprian  and  Tu  (1983)  and  Mizushina, 
Maruyama  and  Hirasawa  (1975)  in  periodic  pipe  flow.  Detailed  studies  showed  that 
the  phase  lag  of  (u3)  (and  (—uv)  )  increases  approximately  linearly  through  the  un¬ 
steady  vortical  layer.  The  flow  thus  behaves  as  if  the  disturbance,  imposed  by  the 
oscillation  on  the  turbulence  generated  near  the  wall,  diffuses  outward  with  a  nearly 
constant  velocity.  The  distance  travelled  by  the  disturbance  in  one  oscillation  period 
due  to  this  diffusion  corresponds  to  the  thickness  of  the  unsteady  vortical  layer.  If  this 
thickness  can  regarded  to  be  characterized  by  ym  (in  order  of  magnitude)  the  present 
data  suggest  that  the  diffusion  velocity  is  of  the  order  of  u..  Another  significant  point 
to  be  noted  is  that  both  (u3)  and  (— uv)  are  out  of  phase  with  (17)  over  most  part 
of  the  unsteady  vortical  layer,  a  result  which  is  in  agreement  with  the  observations  of 
Ramaprian  and  Tu  (1983).  This  points  to  the  limitations  of  quasisteady  eddy  viscosity 
models  that  relate  (— uv )  to  ( U ).  A  careful  examination  of  figure  7  also  shows  that  (u3) 
and  (— uv)  have  a  relative  phase  difference  which  varies  across  the  unsteady  vortical 
layer.  This  is,  again,  in  agreement  with  the  earlier  results  from  pipe  flow  studies  and 
shows  that  the  turbulence  structure  is  not  in  "equilibrium”  during  the  oscillation  cycle. 

The  mean  value,  as  well  as  the  amplitude  and  phase  of  the  first  harmonic  of  the 
skin-friction  coefficient  (Ct)  as  a  function  of  oscillation  frequency  are  shown  in  figure  8. 


These  data  were  obtained  from  a  set  of  experiments  in  which  the  HFG  was  located  at 
station  5  and  the  frequency  of  oscillation  was  varied  from  0  to  2  Hz  by  using  the  sleeves 
designed  for  0.5  Hz  and  2  Hz.  It  was  found  that  there  was  very  little  distortion  in 
the  wave  form  of  the  freestream  velocity  variation  in  spite  of  this  off-design  operation. 
A  slight  decrease  in  the  time-mean  value  C,  (shown  as  a  shaded  band)  is  seen  as 
the  oscillation  frequency  is  initially  increased  but  the  Ct  values  seem  to  level  off  at 
higher  frequencies.  The  amplitude  of  (Ct)  is  normalized  in  the  same  manner  as  (— uv). 
Hence,  the  difference  between  the  amplitude  and  the  time-mean  value  is  a  measure  of 
the  departure  from  quasisteady  behavior.  It  is  seen  that  the  oscillatory  skin  friction 
shows  only  a  slight  departure  from  quasisteady  response  at  all  frequencies.  However, 
a  monotonic  and  linear  increase  of  phase  lead  of  (Ct)  with  frequency  is  observed.  The 
highest  phase  lead  measured  was  about  38  degrees.  These  phase  lead  values  are  higher 
than  the  value  of  about  10  degrees  observed  by  Ramaprian  and  Tu  (1983)  in  pipe  flow. 
They,  however,  seem  to  support  the  earlier  experimental  results  on  periodic  boundary 
layers  of  Jayaraman,  Parikh  and  Reynolds  (1982)  which  suggest  that  the  phase  lead 
of  (r„)  may  eventually  reach  the  laminar  limit  of  45  degrees  at  very  high  frequencies. 

AN  ASYMPTOTIC  THEORY  FOR  THE  PERIODIC  TURBULENT 
BOUNDARY  LAYER  IN  ZERO  MEAN  PRESSURE  GRADIENT 

The  principal  difficulty  encountered  in  the  analysis  of  turbulent  flows  both  in 
steady  and  unsteady  regimes,  is  the  well  known  closure  problem.  A  combination  of  a 
minimum  of  experimental  information,  dimensional  analysis  and  singular  perturbation 
theory  has  proved  useful,  in  recent  times,  in  providing  quite  general  asymptotic  results 
for  large  Reynolds  numbers,  in  the  form  of  similarity  laws  for  many  categories  of  steady 
turbulent  shear  flows  (see  Tennekes  and  Lumely  1980).  The  strength  of  this  approach 
(and,  perhaps,  also  its  major  limitation)  is  that  it  works  on  the  open  system  of  equa¬ 
tions,  thus  avoiding  the  anomalies  introduced  by  a  specific  turbulence  model.  While  in 
the  case  of  steady  flows  such  analysis  has  merely  reconfirmed  already  well-established 
experimental  knowledge,  it  was  used  in  the  present  work  as  a  predictive  tool  to  obtain 
similarity  laws  for  unsteady  turbulent  flows.  Some  of  these  results  were  then  verified 
using  the  available  experimental  information. 

In  analyzing  the  behavior  of  the  velocity  in  unsteady  (periodic)  flows,  one  can 
search  for  similarity  laws  (e.g.,  law-of-the-wall,  velocity-defect  law,  etc.)  in  the  in¬ 
stantaneous  (ensemble-averaged)  profiles  (e.g.,  Cousteix,  Houdeville  and  Javelle  1981). 
This  approach  has  not  so  far  been  very  successful.  In  the  present  work,  instead,  the 
velocity  profile  is  decomposed  into  its  time-mean  and  oscillatory  components,  and  each 
one  of  them  was  studied  separately.  The  analysis  and  the  results  are  briefly  summarized 
below.  More  details  can  be  found  in  Menendez  and  Ramaprian  (1983,  1986). 

The  analysis  is  started  from  the  ensemble-averaged  unsteady  turbulent  boundary 
layer  equations  for  an  incompressible,  two-dimensional  flow  (Patel  and  Nash  1971) 
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The  particular  form  of  the  periodic  free-stream  velocity  given  by  (l)  is  assumed  with 
both  the  time-mean  value,  U0t  and  the  relative  amplitude  of  oscillation  ,  c  (<  1)  being 
constants.  The  assumption  e  <i  1  allows  a  linear  analysis  in  amplitude  to  be  performed. 
This  means  that  only  motions  at  the  forcing  frequency  are  considered  to  be  significant, 
while  higher  harmonic  oscillations  can  be  neglected.  The  advantage  of  the  complex 
notation  can,  then  be  utilized,  defining 


and  restating  (1)  as 
Note  that  Ult  Vx  and  rt  are  all,  in  general,  complex  quantities. 
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In  a  steady  turbulent  boundary  layer  at  large  Reynolds  number,  there  exist 
two  distinctly  different  length  scales  of  motion  in  the  boundary  layer,  resulting  in  a 
double-layer  structure;  an  outer  layer  characterized  by  the  (local)  mean  boundary  layer 
thickness,  6  and  an  inner  layer  with  length  scale  i//u.  (Milikan  1938).  Furthermore,  as 
the  boundary  layer  thickness,  6,  is  an  ill-defined  quantity,  a  thickness  defined  by 
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will  be  used  instead,  where  9  is  the  momentum  thickness.  Note  that  the  ratio  of  the 
two  length  scales  defines  the  following  Reynolds  number: 


(12) 


The  velocity  scale  for  both  the  layers  is  u,.  Hence  the  time  scales  are  t0  =  A/u,  and 
tt  =  i//uj  for  the  outer  and  inner  layers,  respectively  (Tennekes  and  Lumely  1980). 


In  an  unsteady  turbulent  boundary  layer,  the  imposed  oscillation  introduces  a 
third  time  scale,  fa/-1,  into  the  problem.  Its  relation  to  the  outer  and  inner  time  scales 
can  be  described  by  the  following  two  frequency  parameters 
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In  addition,  a  third  length  scale,  namely  u./w,  appears,  being  associated  with  the 
imposed  oscillation.  This  defines  an  "unsteady  vortical  layer”  whose  thickness  bears 
to  the  outer  and  inner  thicknesses  the  same  relations  as  the  corresponding  time  scales, 
and  is  thus  described  by  the  same  parameters,  i.e.  to  and  a,  defined  above. 

Various  flow  regimes  can  be  distinguished,  according  to  the  values  of  the  fore¬ 
going  parameters.  For  to  <  1,  the  time  scale  to~l  >  t0  and  hence,  significant  unsteady 
vorticity  must  be  expected  in  the  outer  layer.  In  fact,  the  unsteady  vortical  layer  is 
merged  with  the  outer  layer.  On  the  contrary,  since  a  =  R,/to  <  R,  »  1,  i.e.,  u>-1  »  t{, 
the  inner  layer  responds  instantly  to  the  excitation,  i.e.,  in  a  quasisteady  manner.  This 
will  be  called  the  "low  frequency  regime”.  Note  that,  in  particular  case  to  c  1,  the 
outer  layer  will  also  behave  in  a  quasisteady  manner. 

The  other  extreme  situation  arises  when  to  >  R,,  i.e.,  a  <  1,  in  which  case 
to~x  <  tt.  Now,  the  unsteady  vortical  layer  will  be  actually  submerged  in  the  inner 
(viscous)  layer.  The  effects  of  unsteadiness  are,  therefore,  significant  in  the  inner  layer. 
On  the  other  hand,  since  to~  1  ^  *o»  there  is  not  enough  time  for  the  turbulent  flow  in 
the  outer  layer  to  respond  to  the  imposed  oscillation  and  hence  this  layer  will  move  as 
a  rigid  body,  i.e.,  it  will  perform  "  irrotational"  or  "slug-flow"  oscillations.  In  fact,  even 
the  outer  part  of  the  inner  layer  may  exhibit  the  slug-flow  behavior.  This  situation  will 
be  named  "very  high  frequency  regime”. 

The  situation  1  <  ii  <  U,  which  lies  between  the  two  limiting  cases  described 
above,  needs  a  more  careful  study.  Since,  under  this  condition,  the  inner  and  outer  time 
scales  are  such  that  C  w~l  <  tQ,  one  should  expect  that  there  will  be  quasisteady 
flow  in  the  inner  layer,  while  a  slug-flow  is  present  in  the  outer  layer.  Furthermore,  the 
fact  that  1//11,  <  i,  <  A  in  this  flow,  suggests  that  the  unsteady  layer  lies  between 
the  outer  and  inner  layers.  Only  in  this  unsteady  layer  are  the  unsteady  vortical 
effects  significant.  In  fact,  it  has  been  found  useful  in  the  analysis  to  divide  the  regime 
liuCil,  into  two  distinct  subregimes:  an  "intermediate-frequency  regime”,  for 
which  1  <L  to  <  J2i/2  (or  equivalently  Rl/3  <  a  <1  R.)  and  a  "high  frequency  regime”, 
with  Rl/2  C  to  <  R,  (or  equivalently  l<a<  J2i/a).  Note  that  to  =  Rl/a  corresponds 
to  the  case  u/1/a  =  (M<)^a*  i.e.,  the  geometric  mean  of  the  outer  and  inner  scales,  for 
which  to  =  a. 

Each  of  the  above  four  frequency  regimes  has  been  analyzed  by  considering  the 
appropriately  normalized  versions  of  equations  (4-6)  and  using  a  combination  of  order 
of  magnitude  analysis  and  matched  asymptotic  expansions.  This  procedure  is  detailed 
in  Menendez  and  Ramaprian  (1983,  1986).  The  procedure  results  in  the  identification 
of  asymptotic  similarity  laws  that  describe  the  variation  of  the  oscillatory  velocity  in 
the  different  layers  in  the  different  frequency  regimes. 

For  a  free-stream  velocity  given  by  (1),  the  longitudinal  velocity  component  can 
be  written  as 

U  =  V  +  t  (Utt  sinu/t  —  t/ia  cos  tot)  (15) 
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where  Utl  and  U13  are,  respectively,  the  in-phase  and  out-of-phase  components  of  the 
oscillatory  motion.  The  similarity  laws  for  Uxl  and  Uia  resulting  from  the  above  analysis 
are  summarized  in  Table  1.  They  are  valid  in  the  asymptotic  limit  for  R ,  — *  oc.  Hence 
they  should  be  referred  to  as  "weak”  similarity  laws  (Yajnik  1070),  in  the  sense  that  they 
hold  only  to  the  lowest  order  asymptotic  analysis.  Deviations  from  these  asymptotic 
laws  should  be  expected,  in  practice,  at  finite  values  of  J2,. 

The  concept  developed  for  zero  time-mean  pressure  gradient  flows  admits  a 
relatively  simple  extension  to  other  types  of  flow.  An  extension  of  the  present  theory 
to  boundary  layers  in  nonzero  time-mean  pressure  gradient,  as  well  as  to  fully  developed 
flow  in  channels  and  pipes  is  also  described  in  Menendez  and  Ramaprian  (1983). 

Experimental  Results  in  the  Light  of  the  Theory 

The  above  theory  has  been  applied  to  the  present  as  well  as  previously  reported 
experimental  results  of  periodic  turbulent  shear  flows.  Detailed  comparisons  are  re¬ 
ported  in  Menendez  and  Ramaprian  (1983,  1986).  Only  a  few  typical  comparisons  are 
shown  in  figures.  Table  2  shows  the  most  relevant  parameters  for  the  present  data  at 
each  station  and  for  the  two  oscillating  frequencies  studied.  It  is  observed  that  the 
majority  of  the  situations  correspond  to  the  intermediate-frequency  regime.  The  re¬ 
maining  ones  are  in  the  transition  range  between  the  low  and  intermediate-frequency 
regimes.  Therefore,  they  are  expected  to  exhibit,  to  some  extend,  the  main  charac¬ 
teristics  of  both  regimes.  Table  3  shows  a  list  of  the  other  experiments,  together  with 
their  most  relevant  parameters. 

The  cases  in  the  low  or  low/intermediate-frequency  regime  must  be  expected 
to  exhibit  the  "logarithmic  law”,  (see  Table  1).  Figure  9(a)  presents  the  in-phase 
oscillatory  component,  Ullt  in  the  inner  coordinates,  for  experiments  Ml,  M5,  CA,  Tl 
and  Jl.  Data  corresponding  to  intermediate  and  intermediate/high-frequency  regimes 
(experiments  32,  J3  and  J4)  are  also  presented  in  this  figure.  Also  shown  is  the  line 
corresponding  to  the  universal  logarithmic  law  for  the  mean  velocity  profile  in  steady 
flow.  Its  slope  is  used  as  reference.  It  is  readily  observed  that,  there  is  a  range  in  which 
the  in-phase  component  varies  logarithmically  and  with  the  same  universal  slope,  as 
predicted  by  the  theory.  Towards  the  outer  edge,  a  "wake  component”  can  be  identified 
as  for  the  steady-flow  mean  velocity  profile,  but  extending  deeper  into  the  boundary 
layer.  Note  that  this  wake  component  is  very  strong  for  experiment  Jl  due  to  the  time- 
mean  adverse  pressure  gradient.  This  reduces  considerably  the  extent  of  the  logarithmic 
region  (situated  around  y+  2  30).  The  corresponding  out-of  phase  component  Uia,  is 
shown  in  figure  9(b).  It  is  seen  that,  in  fact,  there  is  a  region  of  constant  U13  for  the  low 
and  low-intermediate  frequency  experiments,  in  agreement  with  Table  1.  This  value 
generally  corresponds  to  the  minimum  of  each  distribution,  except  for  experiment  Jl 
for  which  the  region  of  constant  value  seems  to  behave  like  a  "saddle  point".  These 
results  show  that,  like  the  time-mean  flow,  the  oscillatory  motion  in  the  outer  layer  for 
the  low-frequency  regime  depends  strongly  on  the  particular  flow  conditions. 

Table  1  shows  that  in  the  intermediate-frequency  regime,  there  must  exist  an 
outer  region  of  "slug"  flow.  That  this  is,  in  fact,  the  case  has  been  confirmed  by  the 


present  experimental  data  as  well  as  the  data  from  other  experiments.  Immediately 
next  to  this  layer,  and  closer  to  the  wall,  an  intermediate  unsteady  vortical  layer  devel¬ 
ops  where  a  "  velocity-defect  law”  (see  also  Table  1)  must  be  satisfied.  Figures  10(a)  and 
10(b)  show  the  in-phase  and  out-of-phase  velocity  components,  respectively,  for  all  the 
present  experimental  data,  plotted  in  the  coordinates  corresponding  to  the  unsteady 
vortical  layer.  The  error  in  the  asymptotic  results  is  of  order  (Ini?,)-1,  which  corre¬ 
sponds  approximately  to  +  0.2  (see  Table  2).  It  is  therefore  concluded,  that  the  collapse 
is  reasonably  good  for  the  truly  intermediate-frequency  cases,  in  agreement  with  the 
theory.  Deviations  occur  for  the  transitional  low/ intermediate-frequency  cases,  spe¬ 
cially  for  the  out-of-phase  component  and  closer  to  the  outer  edge  of  unsteady  layer. 
In  these  cases,  the  slug-flow  region  has  nearly  disappeared,  so  the  unsteady  layer  ex¬ 
tends  practically  up  to  the  edge  of  the  boundary  layer.  Hence,  the  deviations  must 
result  from  the  interaction  with  the  free-stream.  In  all  the  cases,  the  scatter  for  the 
out-of-phase  component  is  much  larger  than  for  the  in-phase  component.  This  is,  as 
already  mentioned,  due  to  larger  experimental  errors  in  the  measurement  of  U12-  It  is 
significant  to  note  that  the  location  of  the  peak  in  the  in-phase  component,  which  is  a 
measure  of  the  extent  of  the  unsteady  vortical  layer,  corresponds  to  y  2  1. 

The  above  discussion  permits  one  to  construct  a  general,  though  still  speculative, 
picture  of  how  the  oscillatory  components  behave  when  the  frequency  parameter  u 
increases,  as  shown  in  figure  11.  This  figure  refers  to  the  case  of  a  zero  time-mean 
pressure  gradient  boundary  layer.  It  is  observed  that  there  is  a  continuous  transition 
between  the  "log  +  wake”  profile  for  the  low-frequency  regime  and  the  Stokes  (viscous) 
solution  for  the  very  high-frequency  regime. 

Generalization  of  the  Theory  to  Other  Flow  Properties 

Though  the  theory  has  been  developed  for  the  oscillatory  velocity  components, 
the  idea  can  be  extended  to  analyze  other  flow  properties  also.  If  the  in-phase  and 
out-of-phase  components  of  the  turbulence  properties,  namely  (ua),  (va),  and  (uv),  are 
defined  with  complete  analogy  with  those  for  the  velocity,  they  must  scale  with  uj.  To 
be  consistent  with  earlier  definitions  they  are  normalized  with  2 u2.  Their  distributions 
in  the  unsteady-layer  coordinates,  for  the  present  experiments,  are  shown  in  figures 
12(a),  (b),  (c).  With  the  same  admitted  error  as  for  the  velocity  (+  0.2),  the  data  for  the 
in-phase  component  collapse  fairly  well  for  all  the  quantities,  except  for  (ua)  in  the  outer 
part  of  the  unsteady  layer.  Deviations  occur  there  for  the  transitional  low/ intermediate- 
frequency  cases.  As  in  the  case  of  velocity,  this  must  be  attributed  to  interactions  with 
the  free  stream,  including  some  feedback  of  free-stream  turbulence  into  the  boundary 
layer.  Less  satisfactory  is  the  performance  of  the  out-of-phase  components,  especially 
in  the  near-wall  region.  This  can  be  due  to  higher  sensitivity  of  the  turbulent  quantities 
to  the  value  of  a.  The  agreement  improves  when  moving  outwards,  except  for  (ua),  for 
which  the  same  remarks  apply  as  for  the  in-phase  component.  Note  that  the  best  data 
collapse  is  obtained  for  (uv). 

As  in  the  case  of  the  turbulence  properties,  the  in-phase  and  out-of-phase  compo¬ 
nents  of  the  wall  shear  stress  (r9)  must  scale  with  ua.  Figure  13  presents  the  amplitude 
phase  of  oscillation  of  the  wall  shear  stress  as  a  function  of  the  frequency  parameter  <Z>. 


Included  are  the  measurements  at  all  stations  and  at  all  frequencies.  The  collapse  is 
satisfactory  within  the  experimental  error  estimated  to  be  +  0.08.  It  is  observed  that 
the  amplitude  does  not  vary  drastically  with  Q.  There  possibly  exists  a  mild  minimum 
around  w  —  30  to  35.  The  variation  of  the  phase  of  (r„)  with  the  frequency  u>  ex¬ 
hibits  an  approximately  linear  trend,  though  the  distribution  should  eventually  reach 
the  45-degree  line  for  large  w. 

THEORY  OF  THE  USE  OF  SKIN  FRICTION  GAGE 

The  conventional  calibration  procedure  for  the  flush-mounted  hot  film,  when 
used  as  a  skin-friction  gauge,  fails  in  high  frequency  periodic  flows.  A  theory  for  the 
use  of  this  gauge  in  unsteady  flows  has  been  developed.  This  theory  leads  to  the  general 
calibration  formula 

r.  =  +  B)‘  +  +  oAijf  (!•) 

which  replaces  the  conventional  calibration  formula: 

Tm*  =  AE 3  +  B,  (17) 

where  E  is  the  anemometer  output  voltage,  A  and  B  are  the  conventional  calibration 
constants  and  et,  c3  are  constants  that  can  be  obtained  by  a  special  calibration  pro¬ 
cedure.  The  usefulness  and  limitations  of  the  proposed  formula  for  laminar  as  well 
as  turbulent  flows  have  been  established  by  comparing  its  accuracy  against  exact  nu¬ 
merical  solutions  of  the  unsteady  thermal  and  hydrodynamic  boundary  layers  over  the 
film.  The  work  is  described  in  detail  in  Menendez  and  Ramaprian  (1084,  1985).  The 
present  theory  has  been  developed  under  certain  restrictive  assumptions,  described  in 
the  above  references.  While  these  assumptions  are  generally  valid  for  most  of  the  ap¬ 
plications  of  the  skin-friction  gauge,  it  is  necessary  to  verify  that  these  conditions  are 
satisfied  before  applying  the  theory  to  any  specific  problem. 

A  procedure  for  the  calibration  of  the  probe  has  also  been  developed.  Also  a 
simplified  linearized  version  of  the  formula  (16)  has  been  suggested.  This  extends  the 
use  of  the  skin-friction  gauge  to  the  measurements  of  the  energy  spectrum  of  wall- 
shear-stress  fluctuations  in  steady  turbulent  flows,  in  which  such  fluctuations  are  small 
relative  to  the  mean  value. 

While  the  dynamic  effects  of  heat  transfer  to  the  fluid  have  been  fully  taken 
into  account  in  the  present  theory,  a  limitation  of  the  present  theory  is  that  dynamic 
effects  on  the  substrate  heat  transfer  have  been  ignored,  albeit  with  some  justification. 
This  aspect  needs  further  theoretical  study.  It  is  also  very  important  to  conduct  ex¬ 
periments  that  can  lead  to  an  assessment  of  the  accuracy  of  the  calibration  procedure 
recommended  in  this  paper.  For  example,  this  procedure  can  be  tested  by  using  it  for 
measuring  the  instantaneous  wall  shear  stress  in  a  laminar  periodic  boundary  layer,  for 
which  the  exact  solution  can  be  obtained  analytically  or  numerically.  Unfortunately, 
because  of  experimental  limitations,  such  test  could  not  be  conducted  by  the  authors. 


It  is  also  possible  to  verify  the  accuracy  of  the  present  procedure  if  alternative  inde¬ 
pendent  techniques  for  measuring  wall  shear  stress  in  unsteady  turbulent  flows  become 
available  in  the  future. 

COMPUTATIONAL  STUDIES 

The  computational  studies  performed  during  the  period  of  the  present  contract 
included: 

1.  the  development  of  "wall  functions”  (see  Patankar  and  Spalding  (1968) 

2.  the  development  of  a  calculation  procedure  for  unsteady  turbulent  boundary 

layers,  based  on  a  2-equation  model  of  turbulence. 

The  wall  functions  appropriate  for  a  Jit  —  e  type  model  were  obtained  for  unsteady 
boundary  layers  from  the  experimental  data  and  the  theoretical  frame  work  described 
in  the  previous  sections.  These  functions  were  obtained  in  the  form  of  algebraic  expres¬ 
sions  that  fitted  the  near-wall  data  in  figures  10,11  and  12.  Details  of  this  procedure 
and  the  recommended  wall  functions  are  described  in  Menendez  and  Ramaprian  (1984). 

A  finite-difference  calculation  procedure  based  on  the  well  known  k  —  e  model  of 
turbulence  was  developed  for  the  prediction  of  unsteady  turbulent  boundary  layers.  The 
method  is  an  extension  of  the  one-equation  procedure  developed  earlier  by  Menendez 
and  Ramaprian  (1981,  1984)  and  has  been  adapted  for  use  with  wall  functions.  The 
procedure  was  used  for  the  prediction  of  the  present  experiments  in  zero  pressure 
gradient.  Some  typical  results  are  shown  in  figures  14  and  15.  These  refer  to  the 
amplitude  and  phase  of  the  oscillatory  velocity,  and  the  in-phase  and  out-of-phase 
components  of  the  Reynolds  shear  stress.  Results  for  both  oscillation  frequencies  of  0.5 
Hz  and  2  Hz  are  shown.  The  results  shown  correspond  to  stations  2,3,4  and  5  and  are 
compared  with  the  experimental  data.  The  initial  conditions  for  the  calculation  were 
matched  with  the  experimental  data  at  station  1.  Only  the  finite-difference  solutions 
are  shown  and  hence  these  do  not  extend  to  the  wall.  The  distributions  near  the  wall 
can  be  obtained  from  the  wall  functions. 

The  agreement  between  the  predictions  and  measurement  can,  at  best,  be  de¬ 
scribed  as  moderate.  The  failure  of  the  calculation  procedure  to  give  better  predictions 
is  most  likely  due  to  the  inadequacy  of  the  turbulence  model  to  describe  the  turbulent 
flow  at  intermediate  to  high  frequencies  of  oscillation.  Inadequacy  of  the  wall  functions 
and  numerical  deficiencies  may  also  be  partly  responsible  but  are  unlikely  to  be  the 
nuyor  factors.  More  effort  is  obviously  needed  to  improve  our  capability  to  predict 
unsteady  turbulent  boundary  layers. 

EXPERIMENTS  IN  ADVERSE  PRESSURE  GRADIENT 

These  experiments  were  performed  to  obtain  a  comprehensive  data  base  on 
periodic  turbulent  boundary  layers  characterized  by  a  significant  degree  of  flow  rever¬ 
sal  during  the  oscillation  cycle.  These  data  compliment  those  obtained  by  Cousteix, 


Houdeville  and  Javelle  (1981),  Simpson,  Shivaprasad  and  Chew  (1983),  Simpson  and 
Shivaprasad  (1983),  and  Jayaraman,  Parikh  and  Reynolds  (1982).  The  present  exper¬ 
iments  provide,  in  particular,  data  on  the  Reynolds  shear  stress  obtained  from  two- 
component  LDA  and  wall  shear  stress  obtained  from  using  a  heat-flux  gage  (corrected 
for  unsteady  effects)  at  intermediate  to  large  values  of  reduced  frequency.  The  am¬ 
plitude  and  phase  of  the  freestream  velocity  were  maintained  approximately  constant 
over  a  significant  distance  along  the  flow.  Two  sets  of  experiments  were  performed.  In 
one  of  these,  a  "mild”  adverse  pressure  gradient  was  set  up  such  that  an  equilibrium 
boundary  layer  would  be  produced  in  steady  flow.  In  the  other  experiment,  a  stronger 
adverse  press  lire  gradient  was  maintained  over  a  short  downstream  portion  of  the  flow. 
It  was  the  intention  to  produce  zero  mean  skin  friction  in  this  flow.  Unfortunately, 
this  could  not  be  achieved  due  to  apparatus  size  limitations.  It  was,  however,  possi¬ 
ble  to  produce  flow  reversal  during  a  substantial  part  of  the  oscillation  cycle,  in  this 
experiment. 


Modification  of  the  Test  Section 


The  desired  pressure  gradient  was  set  up  in  each  experiment  by  appropriately 
modifying  the  shape  of  the  top  wall  of  the  test  section.  In  the  first  case,  this  shape  was 
calculated  by  imposing  the  condition  that  the  Clauser  parameter  0  defined  as: 

p  ul  dx 

should  remain  constant  in  steady  flow  in  the  downstream  part  of  the  test  section.  The 
actual  constraint  imposed  on  0  is  shown  in  figure  16.  This  distribution  of  0,  along 
with  the  appropriate  initial  conditions,  was  used  in  an  inverse  procedure  to  calculate 
the  required  top  wall  geometry  for  obtaining  an  equilibrium  boundary  layer  in  steady 
flow.  Head’s  entrainment  method  was  used  for  the  boundary-layer  calculation  in  this 
procedure.  The  resulting  top-wall  geometry  and  the  measurement  locations  are  shown 
in  figure  17(a).  The  measured  values  of  the  boundary-layer  parameters  H  and  II  are 
shown  in  table  4.  It  is  seen  that  the  steady  boundary  layer  is  nearly  in  equilibrium 
over  the  region  x=1.42  to  2  meters. 

In  the  second  experiment,  the  shape  of  the  top  wall  was  further  modified  in 
the  downstream  section  in  order  to  obtain  a  higher  pressure  gradient.  The  upstream 
part  of  the  test  section  remained  the  same  as  in  the  first  experiment,  but  a  50  cm 
long  flat  plate  set  back  at  an  angle  of  incidence  of  12  degrees  was  used  to  replace  the 
downstream  part  of  the  upper  wall,  as  shown  in  figure  17(b).  The  boundary  layers  over 
a  length  of  5  cm  along  the  top  and  side  walls  just  upstream  of  the  plate  was  withdrawn 
through  perforations  in  the  wall  to  prevent  flow-separation  along  these  walls.  The 
arrangement  for  the  withdrawal  of  the  boundary  layers  is  shown  in  figure  18.  The 
dimensions  of  the  various  components  of  this  system  were  calculated  from  numerically 
modelled  unsteady  flow  through  the  system,  driven  by  the  oscillating  pressure  inside 


the  tunnel.The  numerical  calculations  indicated  that  the  (two-dimensional)  boundary 
layer  would  separate  inside  the  tunnel  at  all  flow  rates  through  the  tunnel.  Actual 
measurements,  however,  showed  that  flow  separation  occurred  only  at  very  low  flow 
rates  in  steady  flow.  In  the  unsteady  flow,  flow  reversal  was  found  to  occur  only  during 
a  portion  of  the  oscillation  cycle.  Stronger  pressure  gradients  could  not  be  set  up 
because  of  the  limited  size  of  the  tunnel. 


Experimental  Conditions 

All  the  initial  experiments  necessary  for  documenting  the  flow  quality  were  per¬ 
formed  in  exactly  the  same  manner  as  in  the  zero-pressure-gradient  (ZPG)  experiments. 
After  these  initial  testa,  regular,  experiments  were  performed  at  two  oscillation  frequen¬ 
cies,  namely  0.5  Hz  and  2  Hz  with  the  mild  pressure  gradient  and  at  only  one  frequency 
of  0.5  Hz  with  the  stronger  pressure  gradient.  The  average  amplitude  of  oscillation  was 
about  0.4  of  the  mean  velocity  in  all  the  cases.  The  variations  of  amplitude  and  phase 
of  the  free  stream  velocity  along  the  entire  test  section  are  shown  for  the  mild-pressure- 
gradient  experiment  in  figures  19(a)-(d).  It  is  seen  that  while  these  variations  are  not 
large(<  5  %  in  amplitude  and  within  +  2.5  %  in  phase)  over  the  downstream  half 
of  the  test  section,  there  is  a  clear  indication  of  spatial  periodicity  in  the  variations. 
This  is  due  to  the  presence  of  such  a  spatial  periodicity  in  the  displacement  thirlrmHM 
(as  observed  by  Cousteix  et  al  1981),  which  is  transferred  to  the  free  stream  by  the 
interaction  between  boundary  layer  and  the  core  flow. 


Summary  of  Results 

As  in  the  case  of  the  ZPG  results,  all  the  experimental  data  have  been  archived 
on  tape  and  is  available  to  any  interested  user.  Only  a  few  typical  results  are  presented 
here.  A  more  elaborate  presentation  and  discussion  of  the  data  can  be  found  in  the 
forthcoming  Ph.D.  thesis  of  B.  Gajdeczko.  The  data  presented  here  correspond  to  the 
mean  and  fundamental  Fourier  component  only.  Higher  harmonics  (information  on 
which  is  available  on  tape)  were  not  found  to  be  significant  in  the  case  of  velocity. 
While  these  are  significant  in  the  case  of  wall  shear  and  turbulent  stresses,  the  fun¬ 
damental  component  is  still  the  moet  dominant  one  even  in  these  cases.  Results  for 
mild  pressure  gradient  are  summarized  first.  These  are  followed  by  the  results  for  the 
stronger  pressure  gradient. 


Results  for  Mild  Pressure  Gradient 
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Time-mean  flow 

The  time-mean  velocity  distributions  in  steady  and  unsteady  flows  are  shown 
in  figures  20(a), (b),(c).  It  is  seen  that  at  this  equilibrium  mild  pressure  gradient,  the 
time-mean  velocity  is  not  affected  by  the  imposed  unsteadiness.  Effects  due  to  the 
shift  in  the  virtual  origin,  observed  in  ZPG  flow  are  absent  here,  apparently  because 
the  rapid  changes  in  the  tunnel  cross  section  cause  the  transition  to  be  fixed  at  the 
same  location  in  both  steady  and  unsteady  flows. 

The  data  on  long-time  averaged  turbulent  shear  stress  are  shown  in  figures  21(a)- 
(d).  This  is  seen  to  be  affected  by  oscillation,  especially  in  the  downstream  locations  and 
the  higher  oscillation  frequency.  Similar  results  were  obtained  for  the  other  Reynolds 
stress  components  also. 


Oscillatory  flow 

The  amplitude  and  phase  of  the  oscillatory  velocity(fundamental  component) 
are  shown  in  figures  22  and  23.  These  properties  exhibit  very  complex  distributions 
across  the  boundary  layer.  The  most  noticeable  feature  of  these  distributions  is  the 
appearance  of  a  two-layer  structure,  especially  in  the  more  downstream  locations  and 
a  lack  of  self-similarity.  The  amplitude  variations  are  small  at  the  higher  frequency,  as 
expected.  The  lack  of  similarity  is  perhaps  due  to  the  existence  of  a  spatial  periodicity 
in  the  direction.  Such  periodicity  could  not  be  observed  clearly  in  the  present  boundary- 
layer  data,  mainly  because  of  the  sparse  spacing  of  the  measuring  stations  along  the 
tunnel.  The  phase  distributions  also  indicate  a  two-layer  structure,  with  the  outer  layer 
exhibiting  significant  phase  lags. 


The  skin-friction  variation  during  the  oscillation  cycle  is  shown  in  figures  24 
(a)-(d).  These  data  were  obtained  from  the  skin-friction  gage,  using  the  procedure 
described  earlier.  While  the  data  around  ( C t)  «  0  are  suspect,  it  is  expected  that  the 
results  for  (Cf)  >  0.005  are  of  acceptable  accuracy.  Furthermore,  these  figures  show 
that  the  present  procedure  allows  negative  values  of  ( Ct )  to  be  measured. 


Finally,  figures  25  and  26  show  the  distributions  of  the  amplitude  and  phase  of 
(uv)  across  the  boundary  layer.  Comparison  of  the  amplitude  distributions  with  the 
distributions  of  JTB  in  figure  21  shows  that,  except  at  station  2  for  0.5  Hz,  the  response 
of  the  turbulent  stress  significantly  departs  from  a  quasisteady  behavior.  The  smaller 
modulations  of  ( uv )  at  higher  values  of  reduced  frequency  and  the  two-layer  structure  of 
the  flow  in  the  downstream  region  of  the  flow  are  quite  apparent  from  the  figures.  The 
phase  distributions  of  {uv)  in  figure  26  show  abrupt  jumps.  These  are  associated  with 
the  transition  between  the  two  layers  of  the  multilayered  structure.  The  phase  jumps 
suggest  that  the  two  layers  are  driven  by  independent  mechanisms.  The  generation  of 
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maximum  turbulent  energy  far  away  from  the  wall  may  be  the  probable  cause  of  the 
double-layer  structure.  However,  no  definitive  explanation  can  be  given  at  this  time. 


Result*  for  Strong  Adverse  Pressure  Gradient 


The  mean  flow  conditions  for  this  experiment  are  shown  in  table  4.  It  can  be  seen 
that  the  data  have  been  obtained  at  close  spacings  in  the  range  x=1.35  to  1.70  meters. 
Over  this  short  region,  the  shape  factor  increases  from  1.48  to  1.70  (compared  to  about 
1.52  in  the  first  experiment).  The  wake  parameter  II  also  is  correspondingly  larger 
in  this  case.  The  unsteady  flow  was  found  to  exhibit  reversal  over  a  greater  fraction 
of  the  oscillation  cycle  in  this  experiment.  As  already  mentioned,  only  one  oscillation 
frequency,  namely  0.5  Hz  was  studied  in  this  case.  Typical  results  are  presented  in 
figures  27-32. 

The  results  are  qualitatively  similar  to  those  at  the  milder  pressure  gradient.  It 
is,  therefore,  not  necessary  to  discuss  them  in  detail.  The  only  important  points  that 
need  to  be  mentioned  are: 


1.  The  time-mean  velocity  profiles  seem  to  indicate  a  slight  effect  of  imposed  os¬ 
cillation  (The  rather  large  effect  seen  at  x=155  cm  cannot  be  explained  except 
as  being  probably  due  to  measurement  errors). 

2.  The  long-time  averaged  Reynolds  shear  stress  shows  a  significant  effect  of  im¬ 
posed  unsteadiness  in  the  downstream  stations 

3.  The  amplitude  of  oscillatory  velocity  exhibits  even  more  complex  distributions 
that  vary  from  station  to  station.  The  phase  variations  are  qualitatively  similar 
to  those  in  mild  pressure  gradient  but  exhibit  stronger  phase  lags  in  the  outer 
layer. 

4.  The  amplitude  and  phase  of  (uv)  exhibit  trends  similar  to  those  in  mild  pressure 
gradient.  In  fact,  there  is  a  stronger  evidence  of  a  two-layered  structure  in  this 
flow. 

CONCLUDING  REMARKS 

The  Tnain  achievements  of  the  research  work  performed  under  this  contract  can 
be  summarized  as  follows. 


1.  An  extensive  set  of  data  have  been  obtained  on  periodic  turbulent  boundary 
layers  in  zero  as  well  as  adverse  pressure  gradient,  from  experiments  performed 
in  a  specially  designed  water  tunnel.  These  data  include  Reynolds  stresses  and 
wall  shear  stress.  The  initial  and  boundary  conditions  are  all  well  documented. 


All  the  data  have  been  archived  on  magnetic  tape  in  a  "standard”  format  and 
are  available  to  the  research  community. 

2.  A  comprehensive  theory  has  been  developed  for  periodic  turbulent  boundary  lay¬ 
ers.  This  theory  leads  to  the  identification  of  several  similarity  laws  for  unsteady 
turbulent  boundary  layers.  These  laws  have,  in  turn,  been  used  to  develop  "wall 
functions”  needed  in  many  boundary-layer  calculation  codes. 

3.  A  theory  has  been  developed  for  the  use  of  flush-mounted  hot-film  gages  to 
measure  skin-friction  in  unsteady  flows. 

4.  The  large- amplitude  periodic  turbulent  boundary  layer  in  time-mean  adverse 
pressure  gradient  exhibits  a  two-layer  structure.  This  structure  is  not  well  under¬ 
stood  at  this  time.  Work  in  this  direction  will  be  continued  as  a  part  of  the  Ph.D. 
dissertation  of  B.  Gajdeczko.  His  thesis  which  is  expected  to  be  completed  by 
about  June  1987  will  also  contain  more  details  on  the  adverse-pressure-gradient 
boundary-  layer  experiments,  than  have  been  presented  in  this  brief  report. 
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Table  1 

Similarity  Laws  for  Periodic  Turbulent  Boundary  Layers 


Freq.  Regime  Low/Int  Low/Int  Low/Int  Low/lnt  Int 
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Table  4 

Mean-Flow  Properties  for  the  Adverse  Pressure  Gradient 
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Mild  Adverse  Pressure  Gradient 
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Figure  2.  Longitudinal  distribution  of  time-mean  global  properties 
of  the  boundary  layer.  A,  steady  flow;Q,  f  *  0.5  Hz;0, 
f  *  2Hz. 


Figure  3.  Mean  velocity  distributions  in  the  boundary  layer  in  the 
inner  and  outer  coordinates.  (a,b):  0.5Hz  (c,d):  2Hz. 

Symbol,  station:  +  ,!;□, 2;  ^,3;  L, 4;  0,5;  - , 

steady  flow;  — .universal  log  law. 
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4. 


distribution  of  u  ,  v  ,  and  -uv  at  Station  5  (a)  0.5Hz  (b)  2Hz  (0 
u^*.  A,  v2«,  -uv}  (unsteady)  {*,  u?;  +,  »,  -uv}  (steady). 


f>  and  <u  >  during  the  oscillation  cycle,  Station 
u.5Hz;  0,  2Hz.  (ii)  <u2>  (at  2Hz):  (a)  r,  =  1.8 
c)  n  =  0.068,  (d)  n  =  0.023. 
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Figure  6  .  Behavior  of  the  oscillatory  component  of  <U>  at  station  5. 
(a)  amplitude  (b)  phase  *,  0.5Hz;d,  2Hz.  Note  [  ] 
denotes  amplitude  and  {  }  denotes  phase. 
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Figure  8.  Frequency  dependence  of  wall  shear  stress  at 
Station  5.  (a)  amplitude  (b)  phase.  □,  using 

sleeve  designed  for  0.5Hz;  0,  usiruj  sleeve  designed 
for  2.0Hz;  ////,  Time-mean  value  C^. 


Oscillatory  velocity  components  in  the  inner  coordinates  (a)  in-phase 
component  (b)  out-of-phase  component.  A,  experiment  Ml ; D ,  M5;  *,  CA 
0,  T1 ;  V,J1;  o,J2;  +  .J3;  x,J4;  -.universal  logarithmic  law. 


Oscillatory  velocity  components  in  the  unsteady  vortical  layer  coordinates 
Present  experimental  results  (a)  in-phase  component,  (b)  out-of-phase 
compomnent.  9,  experiment  Ml ; 0  *  M2;  M3;  A,  M4;  a,  M5;  #,  M6;  X,  M7; 
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Figure  11, 


Qualitative  representation  of  the  evolution 
of  the  in-phase  and  out-of-phase  velocity 
components  with  w.  0<  <  u2  <  . ..<  u7. 

(a)  In-phase  component,  (b)  Out-of-phase 
component. 
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Figure  12,  Oscillatory  components  c 
<u^>  <v^>  and  -<uv>  in 


unsteady  vortical  layer 
coordinates.  Symbol  as  in 
figure  11.  Subscript  11 
denotes  in-phase,  subscript 
22  denotes  out-of-phase. 
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Figure  13.  .Correlation  of  the  oscillatory  wall  shear  stress  in  terms 
of  w.  (a)  Amplitude,  (b)  Phase.  □  ,  f  -  0.5  Hz.O,  f  •  2 
a,  station  5,  0.5  Hz  sleeve;  v,  station  5,  2  Hz  sleeve; 
O,  quasi -steady  value. 
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Fig.17.  Top  wall  geometries  for  mild  and  strong  pressure  gradients.  The  velocity 
measurement  stations  are  indicated  by  the  arrows.  The  locations  of  HFG  are  marked  on 
the  lower  wall.  For  the  strong  adverse  pressure  gradient,  the  section,  where  boundary 
layer  suction  was  introduced  is  shown. 
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Fig.18.  Detail*  of  the  suction  arrangement  for  the  control  of  the  boundary  layer. 


Fif.26.  Distributions  of  the  phase  of  (-uv)  in  mild  adverse  pressure  gradient 


Fig.27.  Time-mean  velocity  distributions  in  strong  adverse  pressure  gradient 
universal  logarithmic  law. 


Fig .28.  Long-time  averaged  Reynolds  shear  stresses  in  strong  adverse  pressure 


Fig.30.  Distribution*  of  the  phase  of  the  oscillatory  velocity  in  strong  adverse 
prenure  gradient. 


